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INTRODUCTION 


PURPOSE. 

The  purpose  of  this  project  was  (1)  to  evaluate  the  possibility  of  utilizing  the 
Schaefer-Doswe 1 1  technique  to  objectively  analyze  a  3-dimensional  wind  field  in 
the  airport  vicinity  using  parameters  transmitted  through  the  Mode  S  data  link, 

(2)  to  determine  the  observation  and  aircraft  density  required  to  resolve  the  main 
features  of  the  atmospheric  wind  flow  in  the  terminal  radar  coverage  area,  and 

(3)  to  test  whether  this  analysis  could  be  performed  in  a  quasi-real  time  frame. 

BACKGROUND. 

In  today's  air  traffic  control  system,  surveillance  and  communication  are  provided 
respectively  by  the  Air  Traffic  Control  Radar  Beacon  System  (ATCRBS)  and  by  very 
high  frequency  (VHP)  voice  radio  between  the  pilot  and  the  controller.  The  pro¬ 
jected  growth  of  air  traffic  has  necessitated  an  increased  automation  of  the  air 
traffic  control  process.  This  need  is  being  met  through  the  development  of  the 
Mode  S  which  combines  an  improved  beacon  interrogator/transponder  system  with 
an  integral  ground-air-ground  data  link  (reference  2).  This  data  link  will  make 
available  such  aircraft  parameters  as  true  airspeed  and  heading,  which  together 
with  groundspeed  and  track  will  permit  computations  of  the  wind  field  in  the  Mode 
S  coverage  area  when  an  adequate  number  of  aircraft  are  presents  The  Mode  S  data 
Link  will  supply  information  to  a  central  ground-based  computer  which  will  imple¬ 
ment  one  of  the  many  available  objective  analysis  techniques  described  in  the 
meteorological  literature. 

The  objective  analysis  techniques  vary  in  complexity  from  simple  interpolation  to 
more  advanced  methods  involving  the  solution  of  any  number  of  partial  differential 
equations  One  of  the  techniques  available  is  that  of  Schaef er-Doswe 1 1 
(reference  i),  which  combines  features  developed  by  previous  investigators  to  try 
to  control  Lhe  growth  of  the  errors  inherent  in  the  data  collection  and  transmittal 
processes  while  trying  to  conserve  the  divergence  and  vorticity  of  the  obser¬ 
vations.  These  errors  are  propagated  through  the  numerical  processing  used  to 
implement  the  objectiv  analysis. 

Information  supplied  by  the  aircraft  by  means  of  the  data  link  will  permit  the 
calculation  of  the  u  and  v  components  of  the  wind  at  the  aircraft  locations.  From 
these  data,  divergences  and  vorticities  are  calculated  and,  together  with  the  u  and 
v  components  of  the  wind  vector,  interpolated  to  the  desired  grid  pattern.  Through 
the  numerical  solution  of  two  Poisson-type  partial  differential  equations,  an 
improved  wind  field  can  be  produced.  This  method  forces  the  wind  field's 
divergence  and  vorticity  to  be  as  close  to  the  predetermined  values  of  the  data  as 
possible  under  the  constraint  that  the  winds  around  the  boundary  will  correspond 
with  those  obtained  by  some  interpolation  technique.  The  u  and  v  components 
determined  can  then,  if  desired,  be  converted  into  the  magnitude  and  direction  of 
the  wind  vector.  The  implementation  of  this  numerical  procedure  will  take  approxi¬ 
mately  two  minutes  of  processor  time  and  30,000  words  of  memory  for  most 
mi nicomputers . 


DISCUSSION 


BRIEF  TECHNIQUE  DESCRIPTION. 

The  numerical  procedure  used  Co  compute  an  objectively  analyzed  wind  field  from 
known  aircraft  observations  proceeds  as  follows: 

].  Calculate  wind  components  at  given  time  intervals  along  the  aircraft  track 
by  means  of  Che  data  supplied  through  the  Mode  S  data  link;  aircraft  positions  and 
wind  components  are  the  only  information  needed  to  implement  the  numerical  method. 

2.  Select  triangles  with  vertices  at  three  aircraft  positions  by  taking  all 
aircraft  positions  successively,  three  locations  at  a  time,  to  form  the  triangles. 
Equilateral  triangles  with  smallest  possible  area  are  preferred  for  the  calculation 
of  che  divergences  and  vorticities. 

3.  Using  a  Bellamy  technique  (reference  3),  or  equivalent,  compute 
divergences  and  vorticities  for  each  of  the  triangles  selected  and  assign  them  to 
the  triangles'  centroids. 

4.  Interpolate  wind  components,  divergences  and  vorticities  to  grid  points 
using  a  convenient  interpolation  technique,  such  as  the  method  described  by  Barnes 
( reference  4) . 

5.  Solve  two  Poisson-type  partial  differential  equations  with  Dirichlet-type 
boundary  conditions.  The  values  of  the  function  on  the  boundary  are  obtained  from 
the  interpolation  technique.  The  Poisson  equations  are  the  Eul er-Lagrange 
equations  obtained  from  a  variational  formulation  (described  later).  The  partial 
differential  equations  were  solved  by  a  finite  difference  procedure,  which  uses  a 
second  order  analogue  for  the  Laplacian  operator  (reference  3),  the  solution  will 
produce  a  correction  on  the  original  interpolated  wind  field. 

Since  Mode  S  and  its  data  link  are  not  operationally  available,  we  have  introduced 
a  known  wind  field  to  compare  with  the  predictions  of  the  mathematical  model.  This 
known  field  was  obtained  by  using  an  anlytical  expression  for  the  u  and  v  com¬ 
ponents  of  the  wind  field,  either  directly  or  through  the  calculation  of  the  field 
components  from  the  stream  function  and  velocity  potentials,  since  the  horizontal 
wind  v  can  be  partitioned  into  irrotational  and  nondivergent  components 

v  =  V[|  *  *  k  x  Vu?  ( i ) 


where  X  is  the  velocity  potential,  f  the  stream  function  and  the  horizontal 

gradient  operator.  The  analytic  expressions  were  selected  so  as  to  model  realistic 
atmospheric  disturbances. 

TECHNICAL  DESCRIPTION  OF  THE  NUMERICAL  METHOD. 

The  problem  involves  minimizing  the  difference  between  the  final  computed  field 
and  the  interpolated  field,  such  that  the  final  field  will  have  the  measured 
divergence  and  vorticity.  One  way  of  doing  this  is  to  cast  the  problem  into 
one  of  a  variational  type  (reference  6),  the  functional  can  then  be  chosen  as 
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J  =  SSs  [(v^-”v^)2  +  A  j  (  k  .  V  X  _V  -  5  )  +  A  2  (V-X-d^  ds  (2) 


where : 


dS  is  the  differential  surface  area,  v^x,y)  final  analyzed  horizontal  vector  wind 
field,  'v(x,y)  preliminary  horizontal  vector  wind,  T(x,y)  measured  vertical  com¬ 
ponent  of  vorticity,  D(x,y)  measured  horizontal  divergence,  and  \  j>A2  Lagrange 
mul t ipl i ers . 


After  setting  the  first  variation  equal  to 

zero  and  using  Green's  Theorem,  we  get: 

v  -  "v"  =  1/2  [V  1  2  +  k  x  ^7  A^  1 

(3) 

)<=> 

ii 

>1 

t> 

(4) 

k  .  V  x  _v  =  i 

(3) 

These  are  the  Euler-Lagrange  equations, 
conditions  are: 

The  corresponding  natural  boundary 

f  '  1  6  v  .  dt  =  0 

(6) 

J  dS  -  - 

fas 

(7) 

where  d  S  is  the  boundary  of  the  surface  and  £  is  the  tangential  unit  vector  along 
the  boundary  curve. 

From  the  Euler-Lagrange  equations,  the  Lagrange  multipliers  satisfy  the  following 
equat i ons : 


V2  A  !  =  2[  1  -  k.  V  x  v] 


(8) 


V  2  X  2  =  2 1 D  -  V-v] 


(9) 


It  is  necessary  to  solve  these  two  partial  differential  equations  with  the  above 
natural  boundary  conditions. 

The  simplest  specification  of  boundary  conditions  that  satisfy  the  natural  boundary 
conditions  would  be: 


X  1 


ds 


(10) 
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Schaefer-Doswell  (reference  1)  shows  Chat  this  technique  gives  a  statistically 
improved  wind  field  analysis,  as  compared  with  ordinary  interpolation  techniques 
and,  at  the  same  time,  preserves  the  original  divergence  and  vorticity. 


This  method  requires  the  solution  of  two  Poisson  equations  using  numerical 
differentiation  to  find  the  grid-point  wind  components.  Although  the  method  is 
fairly  accurate,  it  is  computationally  complicated.  We  will  use  a  second  method 
which  is  computional ly  less  costly;  this  method  produces  a  wind  whose  divergence 
and  vorticity  are  as  close  to  those  measured  as  possible  under  the  restriction  of 
boundary  specification  of  the  winds. 

A  new  variational  principle  is  now  introduced  which  does  not  include  the 
difference  Cv  -  v.)2  ;  i.e.,  the  restriction  that  the  difference  between  preliminary 
and  analyzed  field  be  minimal  is  dropped.  The  functional  is  given  by: 

1  =  //s  l(k.  V  x  v  -  T)2  +  (7.  v  -  D)2|  dS  nn 


After  taking  the  first  variation,  this  becomes: 
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where  m  is  the  metric  coefficient  determined  by  the  analysis  grid  geometry. 
The  natural  boundary  conditions  will  be  satisfied  when  either 


(V.  v  -  D) 


gs  =  (k.Vxv  -  () 


0 


(14) 


or 


(15) 


Because  of  the  interrelationship  between  V.  _v  and  k  .  V  x_v,  the  second  boundary 
condition  would  be  appropriate;  this  condition  will  be  satisfied  if  the  winds  along 
the  boundary  are  specified. 


When  this  technique  is  compared  with  interpolation  techniques,  a  marked  improvement 
is  shown  as  the  domain  of  comparison  is  moved  away  from  the  boundaries 
(reference  7).  Therefore,  the  quality  of  the  wind  analysis  is  enhanced  by 
increasing  the  accuracy  and  number  of  the  wind  observations  along  the  boundary. 
Nevertheless,  when  the  outer  grid  values  are  omitted,  a  marked  improvement  occurs 
with  the  variational  technique  over  the  interpolation  technique. 

DETAILED  DESCRIPT10 N  OF  COMPUTATIONAL  PROCEDURE . 

After  computing  the  wind  velocity  components  at  the  aircraft  locations  along  the 
track,  divergences  and  vorticities  must  be  computed  from  the  wind  velocity  com¬ 
ponents.  The  method  described  by  Schaefer-Doswe 1 1  will  force  the  analyzed  wind 
field's  divergence  and  vorticity  to  be  as  close  to  the  predetermined  values  as 
possible  under  the  restriction  that  winds  around  the  grid  boundary  will  match  those 
obtained  through  interpolation  of  the  original  wind  observations. 

The  horizontal  divergence  D  can  be  computed  from  the  definition: 

D  V  •  v  =  3u  3v  (16) 

H  3x  3y 

This  evaluation  of  D  through  the  differential  definition  would  cause  two  problems: 
(1)  interpolation  of  winds  to  a  uniform  grid  is  a  nonunique  process  and  (2)  a 
centered  differences  computation  would  underestimate  the  derivative.  Morel  and 
Necco  (reference  8)  have  shown  that  the  total  uncertainty  of  the  computation  can 
exceed  100  percent.  To  avoid  some  of  these  problems  (grid  interpolation  would  no 
longer  be  necessary;  discontinuities  inherent  in  interpolation  techniques  would  be 
replaced  by  the  smoothing  function  of  integration),  the  integral  definition  of 
divergence  can  be  used: 


D  1  im  1 
A-0  71 


k . v  x  dr 


A  is  the  area  of  the  surface  S.  We  integrate  around  the  contour  of  the  surface, 
where  3  S  is  the  curve  that  bounds  S.  Bellamy  (reference  3)  has  developed  a  tech¬ 
nique  to  estimate  divergences  using  the  integral  definition.  For  the  contour  3  S, 
we  take  a  triangle  with  vertices  at  the  wind  observation  locations  given  by  the 
aircraft  position  along  the  aircraft's  track.  The  winds  are  then  allowed  to 
displace  the  Locations  of  the  vertices  for  a  small  time  interval  6  t.  The  diver¬ 
gence  then  equals  the  percentage  change  in  area  enclosed  by  the  curve-per-uni t 
time  and  is  assigned  to  the  triangle  centroid.  Similarly,  the  vertical  component 
of  vorticity  £  can  be  defined  by: 


k.  V  xv 


5v  _  lim  _1  I.  v.  dr_ 

) x  dy  "  A*0  A  / 


This  can  be  evaluated  in  Che  same  manner  as  the  divergence,  since  rotating  the  wind 
vectors  by  90°  to  the  right  will  generate  a  wind  field  with  vorticity  equal  to  the 
divergence  of  the  original  field. 


To  construct  the  triangles  required  for  the  calculation  of  the  divergences  and 
vorticities,  a  circle  is  drawn  centered  around  each  observation  sequentially.  The 
radius  is  varied  until  between  five  and  ten  observations  are  included  within  the 
radius  (this  number  is  arbitrary  and  can  be  changed  when  convenient  for  the  calcu¬ 
lations).  All  possible  combinations  of  three  points  at  a  time  are  considered,  and 
the  triangles  are  constructed  from  these  three  observations.  The  triangles  with 
smallest  area  and  closest  to  being  equilateral  are  retained  for  the  calculations; 
i.e.,  the  triangles  are  weighted  acording  to  this  criteria,  and  two  more  triangles 
than  the  number  of  observations  within  the  circle  are  arbitrarily  retained.  After 
che  line  integrals  are  evaluated  and  a  value  is  obtained  for  the  divergence  (or 
vorticity),  this  value  is  assigned  to  the  triangle  centroid.  Values  of  the  diver¬ 
gences  and  vorticities  at  the  triangle  centroids  and  the  values  of  the  wind 
components  at  che  aircraft  observation  locations  are  then  interpolated  to  the  grid 
points  where  the  wind  values  are  desired.  The  grid  spacing  should  ideally  be  of 
dimensions  about  half  the  average  data  spacing;  however,  in  some  situations  the 
data  and  grid  spacing  are  fixed  for  the  particular  problem  and  cannot  be  varied  at 
will. 

The  next  step  consists  of  the  interpolation  of  values  to  grid  points  in  a  regular 
array.  Any  of  the  known  interpolation  procedures  could  be  used  at  this  stage; 
however,  we  employed  a  technique  described  by  Barnes  (reference  4).  This  inter¬ 
polation  has  a  response  function  D(a,k)  =  exp  (-a2  k),  where  a  =  2  tt/  A ,  X  =  wave¬ 
length  and  k  is  an  adjustable  parameter.  Barnes  also  considers  interpolation  in 
time  with  an  exponential  response,  which  we  have  omitted  because  time  weighting  was 
determined  to  introduce  small  errors  when  compared  with  other  sources  of  error. 

When  the  grid  values  are  obtained,  the  partial  differential  equations 


m2  \A,  (' 

i1)] 

-  2  [  d 

-  — 

(Q  ♦  A  (£Y| 

[ dx2  \  n 

o/  3y2 
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Ldx 

\m2/  3y  \m7J 

are  to  be  solved.  (m  is  che  metric  coefficient  determined  by  the  grid  geometry.) 
As  mentioned  before,  these  equations  were  obtained  from  the  variational  functional: 


I  =  Sh  l  ( k .  V  x_v-£)2  +  (V.  v 
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by  means  of  the  Euler- Lagrange  equations.  The  natural  boundary  conditions  are  now 
satisfied  if  either 


(V-v  -  D)  ds  =  (k.Vxv  -  £)  as  =  0 


(22) 


or 


6  *  |  9S  =  0 


(23) 


Because  of  the  interrelationship  between  *7-  v  and  k .Vx  v,  the  second  condition 
will  be  used.  This  condition  can  be  satisfied  if  the  winds  along  the  boundary 
are  prescribed  by  an  independent  technique,  such  as  interpolation. 


To  solve  the  partial  differential  equations,  we  have  used  an  "Improved  Second 
Finite  Difference  Analogue  for  the  Laplacian  Operator,"  Schaefer  (reference  5). 
This  finite  difference  is  given  by 


Af 
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Ffx  +  6,  v+  «S ) 


(x+6,y)  +  F  (x-6,y)  +  F(x,y+6)  +  F(x,y-<5)j 
+  F(x+6,y-i)  +  F(x-6,y-b)  +  F(x-<5,y+6)  - 
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F(x , y) 


(24) 


(6  is  the  grid  spacing).  This  approximation  has  the  property  that  for  a  general 
function  there  is  a  minimum  of  second  order  error  and,  if  the  function  is  defined 
only  at  specified  points,  this  is  the  "best"  possible  approximation  to  the  differ¬ 
ential  operator  (acts  to  filter  out  high  frequency  computational  wave  modes  in  the 
computed  Laplacian  field  by  increasing  the  correlation  between  point-by-point 
values  calculated  by  the  operator  and  the  Laplacian). 

By  solving  these  partial  differential  equations  with  the  boundary  conditions 
obtained  by  the  interpolation  scheme,  an  improved  wind  field  is  produced  with  a 
divergence  and  vorttcity  chat  will  be  as  close  as  possible  to  Che  original 
divergence  and  vorticity,  while  the  winds  around  the  boundary  will  approximate 
those  obtained  by  the  interpolation  technique. 

NUMERICAL  SIMULATION  OF  WIND  FIELDS 

Several  numerical  simulations  were  carried  out  to  test  the  validity  and  accuracy  of 
the  procedure  in  realistic  atmospheric  situations.  We  constructed  analytic  wind 
fields,  sampled  them  at  the  aircraft  positions  around  a  hypothetical  airport,  and 
these  sampled  wind  fields  were  processed  through  the  numerical  scheme  to  be 
compared  with  the  original  analytic  fields. 
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We  tested  analytic  wind  fields  which  consisted  of  one  or  more  waves  of  different 
amplitudes,  periods  and  orientations,  and  also  constructed  analytic  wind  fields 
from  mathematical  expressions  for  stream  functions  that  resemble  atmospheric  flows 
(reference  9);  i.e.,  of  the  form 


V (x , y )  =  C  -  Uy  -  V y2  +  Wy3  + 
-  D  |  1+  (x-x0)2  +  (y-yo)2j/ 


Asinx  +  Bcos2x 


Ecos  x  sinv 


(25) 


where  C,  U,  V,  W,  A,  B,  E,  D,  L,  x  and  y  are  constants.  This  was  used  by 
Schaefer  and  Doswell  to  test  the  divergence  calculations  since  a  field  with  zero 
divergence  should  be  produced  from  a  stream  function. 

We  will  comment  on  and  show  graphically  the  wind  reconstruction  for  a  few 
representative  cases  which  will  illustrate  the  effect  of  changing  some  of  the  more 
influential  parameters;  such  as,  the  effect  of  reducing  the  number  of  observational 
aircraft  within  the  range  of  t he  Mode  S  t rack i ng  antenna,  the  effect  of  changes  in 
sample  rates  per  aircraft,  the  effect  of  observational  errors,  the  effect  of 
atmospheric  field  wavelengths  on  wind  field  reproduction,  etc. 

First,  we  will  consider  a  number  of  aircraft,  (the  number  of  aircraft  will  be 
reduced  from  five  to  two)  collecting  data  at  a  rate  of  one  per  20  seconds  for  a 
period  of  2  minutes  and  flying  through  an  atmospheric  environment  (figure  1) 
where  the  u  component  of  the  wind  (figure  2A)  consists  of  a  mean  value  of  30 
units  with  two  waves,  amplitude  three  units  and  amplitude  one  unit  respectively 
(units  chosen  are  irrelevant  to  the  numerical  procedure),  and  period  of  one  wave 
per  25  grid  points  (grid  spacing  could  be  considered  1  mile)  (figure  2F). 

The  numerical  method  shows  that,  although  the  mean  value  of  the  field  is  regained 
by  the  numerical  reconstructive  procedure,  the  shape  of  the  wave  would  degenerate 
rapidly  as  the  number  of  aircraft  is  reduced  (figure  2B  to  E).  This  rapid  degener¬ 
ation  of  the  wave  shape  is  inherent  in  the  numerical  procedure,  since  to  calculate 
the  divergences  and  vorticities,  a  random  placement  of  the  observations  is  more 
useful  than  the  alignment  of  observations  in  a  reduced  number  of  aircraft  tracks. 
The  shape  of  the  wave  is  no  longer  evident  when  the  number  of  aircraft  falls  below 
three,  if  there  is  a  prominent  disturbance  where  the  field  has  not  been  sampled. 
The  improved  accuracy  of  the  field  as  we  move  away  from  the  edges  is  also  visible 
in  the  figure.  When  the  number  of  aircraft  is  low,  a  straight  interpolation  pro¬ 
cedure  would  be  more  efficient,  since  conservation  of  divergence  and  vorticity 
would  not  appreciably  enhance  the  quality  of  the  reconstruction. 

Next,  we  checked  the  effect  of  increasing  the  data  rate  from  20  seconds  (every 
fifth  antenna  rotation)  to  10  seconds  and  then  5  seconds,  data  rates  that  might  be 
achieved  from  Mode  S  if  required.  We  used  the  same  wind  field  and  considered  the 
effect  of  the  increased  data  rate  with  three  aircraft  collecting  wind  information 
(figure  3).  As  shown,  the  reconstruction  of  the  field  does  not  improve  markedly, 
and  therefore  a  slower  sample  rate  of  one  observation  every  20  seconds  is  suf¬ 
ficient.  The  effect  of  the  number  of  aircraft  is  more  marked  chan  the  number  of 
observations  per  aircraft,  even  when  the  total  number  of  observations  is  the  same, 


8 


nn 

mm 

WmwV'M 

mWM  'M 

wmmm 


'pa 

\u. 

J  W"\l 
Kv  V''l 

r*  *\V  I 

jffMy  •:v'' 

Vrj&M 


mm 

mmM 

|1 


K<wX 

PtViVmi 

Rmmvl 

vtVtVfoVmlv 

mmmm 

P:#M::Svk« 

BpMMfl 

kwx»M 


mm 


mm 

w 


mmm 


ptit|| 

tmm 


SsssMI 


KM 

ftt] 

JBw'wl 
pKvviv>4 
,y$5B 

RW^^AVi 

IM 

ml 

In 

BkwmsSt 

Hr 


Kl 

Pill 

l?#l 

ii 

E&Ml 


EFFECT  OF  REDUCTION  IN  NUMBER  OF  AIRCRAFT  —  SIMULATION  GRID 


ONE  OBSERVATION  PER  20  SECONDS 


ONE  OBSERVATION  PER  10  SECONDS 


ONE  OBSERVATION  PER  5  SECONDS 


FIGURE  3. 


INCREASE  IN  DATA  RATE 


since  the  locations  of  the  observations  are  more  randomly  placed  in  this  case  and 
computations  of  divergence  and  vorticity  become  more  reliable  (there  are  more 
nearly  equilateral  triangles). 

To  show  the  effect  of  errors  in  the  data  on  the  numerical  scheme,  we  introduced 
both  a  random  error  on  all  the  observations  and  a  bias  on  the  observations  from  one 
aircraft.  The  random  errors  are  reduced  about  a  third  due  to  the  smoothing  carried 
by  the  numerical  method.  When  a  bias  of  three  units  on  one  of  five  aircraft 
(figure  4)  is  introduced,  we  see  that  even  in  the  worst  possible  situation  when  the 
aircraft  is  flying  perpendicular  to  the  wavefronts,  only  the  field  in  the  immediate 
vicinity  of  the  biased  airplane  is  appreciably  changed.  Even  this  bias  is  somewhat 
smoothed  at  the  worse  locations  near  the  aircraft  track. 

If  the  number  of  waves  per  grid  area  is  increased,  the  quality  of  the  wind  field 
reproduction  deteriorates  (reference  10).  This  is  a  consequence  of  the  relation¬ 
ship  between  the  mean  observation  separation  and  grid  spacing.  The  grid  spacing 
should  be  chosen  to  be  of  dimensions  roughly  half  that  of  the  average  data  spacing 
with  the  consequent  limitation  on  the  resolvable  frequency  of  the  waves. 


CONCLUSIONS 


The  S c h ae f e r - Do s we  1 1  technique  was  utilized  to  investigate  wind  field 
reconstruction  from  sparse  observations  and  to  simulate  the  effects  on  the 
numerical  technique  of  such  factors  as  increased  data  rate,  changes  in  number  of 
aircraft,  random  and  biased  errors,  and  frequency  of  waves.  In  spite  of  the 
possible  degradation  due  to  other  than  ideal  conditions,  it  seems  that  in  most  cir¬ 
cumstances  this  technique  produces  a  good  visualization  of  the  atmospheric 
environment  and  an  adequate  numerical  reconstruction  of  the  wind  field. 

Other  techniques  are  available  which  can  perhaps  give  better  results,  such  as  the 
development  of  the  first  variational  technique  described;  however,  most  of  the 
techniques  considered  are  costly  in  both  computer  time  and  storage,  while  the 
second  variational  technique  used  seems  to  be  adequate  for  a  general  description  of 
the  airport  wind  environment. 


RECOMMENDATION 


While  our  simulation  work  and  that  of  Schae f er-Doswe 1 1  in  the  synoptic  scale  tend 
to  indicate  that  an  adequate  description  of  the  wind  field  is  possible  through  this 
numerical  technique,  it  would  be  useful  to  provide  actual  measured  wind  fields  to 
the  technique.  As  soon  as  the  Mode  S  data  link  is  available  and  aircraft  can  be 
equipped  with  Mode  S  transponders  and  wind  sensors,  it  would  be  desirable  to  con¬ 
tinue  the  technique  testing.  In  this  way,  the  behavior  of  the  numerical  method 
under  actual  conditions  can  be  tested,  and  corrections  to  the  numerical  method  or 
further  tuning  can  be  implemented  if  necessary.  If  more  accuracy  is  desired, 
perhaps  the  first  variational  technique  can  be  developed  within  our  limitations 
of  processor  time  and  core  requirements. 
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CALCULATED  FIELD  (NO  BIAS  AND  BIAS) 
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